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Abstract 

We consider an initial boundary value problem for a 2x2 system of conservation laws mo- 
deling heatless adsorption of a gaseous mixture with two species and instantaneous exchange 
kinetics, close to the system of Chromatography. In this model the velocity is not constant 
because the sorption effect is taken into account. Exchanging the roles of the x, t variables we 
obtain a strictly hyperbolic system with a zero eigenvalue. Our aim is to construct a solution 
with a velocity which blows up at the corresponding characteristic "hyperbolic boundary" 
{t = 0}. 

AMS Classification: 35L65, 35L67, 35Q35. 

Key words: hyperbolic systems, conservation laws, weak solutions, Temple systems, blow 
up, boundary conditions, gas chromatography, Riemann problem, Front Tracking Algorithm. 

1 Introduction 

This paper deals with the construction of a blow up solution for a one dimensional 2x2 hyper- 
bolic system of conservation laws related to a particular isothermal gas-solid chromatography 
process, called "Pressure Swing Adsorption" (PSA), with two species. Each of the two species 
simultaneously exists under two phases: a gaseous and movable one with a common velocity 
u(t,x) and concentration Ci(t,x) > and a solid (adsorbed) with concentration q^, i = 1, 2. 
One may consult [5D] and [53] for a precise description of the process and [5] for a survey on 
various related models. 

In gas chromatography, velocity variations accompany changes in gas composition, es- 
pecially in the case of high concentration solute: it is known as the sorption effect. Here, 
the sorption effect is taken into account through a constraint on the pressure that is, in this 
isothermal model, on the total density: 

c 1 + c 2 = p(t), (1) 

"Universite de Savoie, LAMA, UMR CNRS 5127, 73376 Le Bourget-du-Lac, bourdarias@univ-savoie.fr 
tT Jniversite de Savoie, LAMA, UMR CNRS 5127, 73376 Le Bourget-du-Lac, gisclon@univ-savoie.fr 
tUniversite de Nice, Labo. JAD, UMR CNRS 6621, Pare Valrose, 06108 Nice, junca@unice.fr 



where the function p is given (this is actually achieved in the experimental or industrial 
device). This constraint expresses that the total pressure depends only on the time. In the 
sequel we assume that p = l. 

We assume that mass exchanges between the mobile and the stationary phases are in- 
finitely fast (instantaneous exchange kinetics) thus the two phases are constantly at com- 
position equilibrium: the concentrations in the solid phase are given by some relations 
q L = q*(ci,c 2 ) where the functions q* are the so-called equilibrium isotherms. A theo- 
retical study of a model with finite exchange kinetics was presented in [5] and a numerical 
approach was developed in [3] . 

With these assumptions, the PSA system reads: 

dt{ci + qt{c u c 2 )) + d x {ucx) = 0, (2) 
dt(c 2 + q* 2 {c u c 2 )) + d x (uc 2 ) = 0, (3) 
ci+c 2 = 1. (4) 

Notice that ci, c 2 must satisfy < c\, c 2 < 1. 
Adding (0) and © we get, thanks to (g]): 

d t (ql(ci,c 2 ) + g|(ci,C2)) + d x u = 0, 

i.e. the constraint (fTJ) leads to an integral dependency of the velocity upon the concen- 
trations. We denote c = c\ then c 2 — 1 — c and the unknowns are the velocity u and the 
concentration c. We write the PSA system under the form: 

d x u + d t h(c) - 0, (5) 
d x (uc) + d t I(c) = 0, (6) 

with 

h(c) = ql(c,l-c)+q2(c,l-c), 
1(c) = c + ql{c,l-c). 

In the sequel we denote 

9i ( c ) = 9i( c > 1 - c) and q 2 (c) = g|(c, 1 - c), 

thus h = qi+q 2 . Any equilibrium isotherm related to a given species is always non decreasing 
with respect to the corresponding concentration and non increasing with respect to the others 
(see [T7]) i-e- 

^f- > and M < for j £ i, 

OCi OCj 

thus we have immediately: 

q'l > > q' 2 - (7) 
PSA system (O-© is completed by initial and boundary data: 

c(0,x) = c (x) e [0,1], .t>0, 

c(t,0) = (%(<)€ [0,1], t>0, (8) 
0) = Ub(t) > a, t > 0, for some constant a > 0. 
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Notice that we assume in ([5} an influx boundary condition, i.e. V< > 0, Ub(t) > and we 
choose ]0, +oo [ instead of ]0, 1[ as spatial domain for the sake of simplicity. 

It is well known that it is possible to analyze the system of Chromatography, and thus 
PSA system, in terms of hyperbolic system of P.D.E. provided to exchange the time and space 
variables: see [21] and also [22] for instance. In this framework, with evolution variable x 
instead of t, the system is strictly hyperbolic and has a null eigenvalue. Thus {t = 0} is a 
characteristic boundary. 

We obtained in [1] [S] an existence result for weak global entropy solutions with large 
BV data for the concentrations and only L°° data for the velocity: as in [T2J [T3] [TH] , the 
zero eigenvalue makes the existence possible of stratified solutions in the sense that u(t, x) — 
Ub(t) v(t, x) where v is as regular as the concentration c and more than the "initial" data it& 
(see [7]). This linearly degenerate eigenvalue makes also possible the propagation of large 
amplitude and high frequency waves as in [111110) . As we will see below, we cannot expect in 
general to have weak global entropy solutions with L°° data for the concentrations because 
it is possible, in that case, to build a blow up solution for some particular isotherms related 
for instance to ammonia or water vapor. 

It is already known that systems of two hyperbolic conservation laws may blow up in the 
L°° norm. In literature ( |27H30j ). there are examples of blow up for one dimensional strictly 
hyperbolic systems of partial differential equations. For 2x2 strictly hyperbolic system 
there are few examples ([IB]). In [3T], see also [35], R. Young built a very nice example 
involving two Burgers equations, linearly coupled at the two boundaries. All examples occur 
in cases where strict hyperbolicity is lost as the solution explodes (see also H.-K.Jenssen 
and R. Young |l6j page 630). In our example, we also loose the strict hyperbolicity of the 
system at the blow up point, but the blow up takes place only at the characteristic boundary 
{t = 0} which becomes twice characteristic and only the velocity blows up. The context of 
our example (gas-solid chromatography) is physically motivated and the various assumptions 
are physically relevant (they are achieved for particular gases) but a blow up along the x-axis 
at t = is of course artificial. Nevertheless it illustrates what may occur when BV regularity 
is not ensured for the velocity at the physical boundary. 

The aim of this paper is to present this example of blow up for a class of PSA systems. 
The paper is organised as follows. In Section [2] we give an overview on PSA system ([S]l-([5| 
and specify our assumptions. In Section [3J we investigate the conditions on PSA system 
©-fS]) to be or not to be in the Temple class: indeed there is no blow up for such systems. 
In Section [4] we state, prove and comment the main result of this paper about a blow up of 
velocity in L°°. 



2 The PSA system 
2.1 General properties 

In this section, we recall briefly some properties of PSA system ([5])-([6j. For more details and 
explanations we refer the reader to [4] [5] [6j [7] . 

As pointed out by Rouchon and al. ([22]), it is possible to analyse PSA system in terms 
of hyperbolic system of P.D.E. provided we exchange the time and space variables. In this 

framework the vector state will be U = \ U ) where to — u c is the flow rate of the first 

V to J 

species. The first component u of this vector must be understood as u p, that is the total flow 
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rate. The reader should bear in mind that convexity is a relevant property of the entropy 
only with respect to the conservative variables (u,m) (see Dafermos [14] section 4.5 page 
76 or [5]). Notice that we seek for solutions satisfying u > 0, wich ensure the system to 
be hyperbolic. Using the diffeomorphism (u,c) i-> (u,m) between R!j_x]0, 1[ and x 
it is more convenient to deal with the hyperbolicity, the Riemann invariants, the Riemann 
problem and the Front tracking algorithm with the set of variables (u, c). For instance there 
is no partial derivative with respect to t for u and c is a Riemann invariant. 
PSA system with conservative variables (u, to) takes the form 



d x U + d t <5>{U) = with $({/) = . (9) 




It is strictly hyperbolic with variable x as the evolution variable: the two eigenvalues are 
(linearly degenerate) and A = H(c)/u where 

H(c) = 1 + q[(c) - ch'(c) = 1 + (1 - c) q[(c) - cq' 2 (c) > 1. (10) 

In particular, the mathematical boundary condition is now related to the set {t = 0} and the 
mathematical initial data to the set {x — 0}. Making use of the function / = q% C2 — q% c\ 
introduced by Douglas and al. in [17], written here under the form 

f(c) = qi {c)-ch(c), (11) 
it is easy to show that the right eigenvector r = ( 1 ^ , \ % ] associated to A satisfies 



l + <7i(c) 

dX ■ r = — ^ /"( c ) 5 thus A is genuinely nonlinear in each domain where /" ^ 0. In view of 
u . 

the remark following (j4]) there is a "natural maximum principle" for the concentration c : 

0<c<l. (12) 

and this is achieved by the solutions (see [1]). 

PSA system ([S])-© admits the two Riemann invariants c and w = ln(u) + 17(c) where g 
satisfies 

Indeed, for smooth solutions: 

d x c + \d f c = 0, 

d x w =0. ^ ' 

In the sequel we make use of the Riemann invariant 

W = uG(c)=e w , 

where 

G(c) = exp( 5 (c)). 
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There are two families of smooth entropies for the PSA system (JS])-©: u ip(c) and <f>(u G(c)) 
where <$> and ^ are any smooth real functions, and the corresponding entropy fluxes satisfy 

Q'(c) =h'{c)^{c)+H(c) i/(c). 

The first family is degenerate convex (in variables (u,m)) provided > 0. So, we seek 
after weak entropy solutions which satisfy, for such functions ip 

0* MM) + t Q(c) <0 

in the distribution sense. 

The second family is not always convex. Nevertheless, if ±G" > then uip = ±uG(c) is 
a degenerate convex entropy, with entropy flux Q = 0, contained in the family of entropies 
uifj(c). So, if G" keeps a constant sign on [0, 1], (c, u) has to satisfy, in the distribution sense: 

±^(uG(c))<0, if ±G" > on [0,1]. (15) 

Notice that the entropies utp(c) and the entropy u G(c) are linear with respect to the velocity 
u. 

2.2 The Riemann Problem 

For the PSA system it is more convenient to begin by the resolution of the boundary Riemann 
problem which is characteristic. The complete Riemann problem is solved at the end of this 
subsection. The boundary Riemann problem has the form: 

d x u + d t h(c) =0, 

d x (uc) + d t i(c) = o, {A>> 
c(0,x) = c°e [0,1], x>o, |^ I £to,' 1] ' * >0 - 

We are classically looking for a selfsimilar solution, i.e.: c(t,x) — C(z), u(t,x) — U(z) with 
z = - > (see Fig. [J). 



A x 




Figure 1: data for the boundary Riemann problem 
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In the domain t > 0, x > 0, the boundary Riemann problem is solved with A— waves 
since A is the only positive eigenvalue of the system. Let us recall the three following results 
obtained in [5], where H and / are given by (fTOl [TTj) . 

Proposition 2.1 ( A rarefaction waves) 

Any smooth non-constant self-similar solution (C{z), U(z)) of H6\) in an open domain 
Q = {0<a<z<f3} where f"(C(z)) does not vanish, satisfies: 



dC_ 

dz 



H{C) 

zf"(cy 



U(z) 



H(C) 



dC 



In particular, —j— has the same sign as /"(C). 

Assume for instance that 0<a<c a <c + <b<l and f" > in ]a,b[. Then the only 
smooth self-similar solution of \lb}) is such that : 



{ C(z) 
dC_ 

dz 
{ C(z) 



H{C) 

Z f"(cy 

c+, 



wh- 



ere z ' = 



U is given by: 



H(c+) 



z a = z+ e -* (c+) with $(c) = / 

Jc° 



< Z < Zq, 

z < z < z + , 
z + < z, 

c no 



(18) 



' U[z) = 
U(z) = 
> U{z) = 



H(C(z)) 



H(0 

< Z < Zq, 



d£. Moreover u° = 



H(c°) 



z 



z < z < z + , 
z + < z. 



(19) 



Proposition 2.2 (A— shock waves) // (c°,c + ) satisfies the following admissibility condi- 
tion equivalent to the Liu entropy- condition ({18}): 

f , 11 W o . + /(^ + ) - f(c°) < f(c) ~ f(c°) 

for all c between c and c , — t < r , 

c+ — c u c — c u 

then the Riemann problem U6\) is solved by a shock wave defined as: 



C(z) 



if < z < s, 
if s < z 



U(z) 



u 



if < z < s, 
if s < z, 



(20) 



where u° and the speed s of the shock are obtained through 



if] 



l + h° 



If] 



l + h+ 



where [c] = c+ - c° , [/] = /+ - f = f(c+) - f(c°), h+ = h(c+), h° = h(c°). 

Proposition 2.3 (A— contact discontinuity ) Two states U° and U + are connected by a 
A— contact discontinuity with c° ^ c + if and only if f is afftne between c° and c + . 
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In this paper, there is no A— contact discontinuity since / is convex: see assumption (H3) 
below. 

It appears from these results that we can build a weak entropy solution of the boundary 
Riemann problem (|16|) in a very simple way (see [5]), similar to the scalar case with flux 
/, for any data. In particular, if /" has a constant sign, the boundary Riemann problem is 
always solved by a simple wave. 

Now, for PSA system (fT6)) . we solve the Riemann problem with the following initial data: 



c(i,0) 
u(t,0) 



c- G [0,1], 



t < 0, 



c(t,Q) 
u(t,0) 



c+ G [0,1], 
u+ > 0, 



t > 0. 



(21) 



u > 0, 

A 0— wave appears on the line {t = 0}. 

Proposition 2.4 (0— contact discontinuity ) Two distinct states U~ and U° are con- 
nected by a 0— contact discontinuity if and only if c~ = c° (with of course u~ ^ u°). 

In conclusion the solution of the Rieman problem for x > and a convex function / is 

• (c, it) = (c~,u~) for t < 0, 

• a 0— contact discontinuity for t = 0, 

• a A— wave for t > 0, 

see Fig. [2] below. In practice, since c° = c~, we first solve the boundary Riemann Problem 
(|16p . (TT71) . Then w° is well defined and the 0— contact discontinuity is automatically solved. 



0-contact disc, 
x 



u 




X- rarefaction 



c , u 



O 



c + , u + 



0-contact disc, 
x 



u 




X- shock 



u 



O 



c + , u + 



c <c + 



c >c + 



Figure 2: solution of the Riemann problem when /" > 0. 



We give also (Fig. [3]) the structure of shock and rarefaction curves when the function / 
is convex, as above, and b! < 0: as stated in [7], these curves are monotonic in that case. 
Moreover the couple of variables (L = Inu, c) is well adapted to the study of these curves 
because L + g{c) is a Riemann invariant. 
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L=ln u 




c 



C" c+ 1 



Figure 3: structure of shock curves (S), rarefactions (R) and contact discontinuities (CD) in the 
case /" > > ti. 

2.3 Main assumptions 

PSA system has the strong property to admit a positive Riemann invariant W = u G(c) > 0, 
which is also an entropy. This Riemann invariant plays a key-role in the blow up mechanism: 
only provided we ensure that W is non decreasing with respect to x we can expect, in view 
of (fl2|) . the velocity u to increase and blow up. Indeed, if W is decreasing with respect to 
x we have < uG{c) < Ub(t) G(cb(t)), thus < u < ||«b||oo supG/ inf G and u is bounded. 

[0,1] P. 1 ] 

Then our first assumption is: 

G" < (HI). 

(HI) means, as seen above (fT5|) . that — W is an admissible degenerate convex entropy for 
weak entropy solutions (with a zero entropy-flux) i.e. d x W > 0. 

In contrary, notice that in some cases, for instance an inert gas and an active gas with a 
Langmuir isotherm as in J4], we have G" > 0, then W is bounded and there is no blow up 
for the velocity in L°° . 

The next hypothesis is useful to the construction of our blow up example in a classical 
hyperbolic framework of Lax with monotonic shock and rarefaction curves as in [7] : 

h' 0, everywhere (H2). 

Assumption (H2) means that one gas is more active than the other. Mathematically, it 
implies the monotonicity of g and G because g' — —h'/H. Thanks to (H2) rarefaction 
curves, which are the level line of W, are monotonic curves in the plane (c, u). 
We also assume that the eigenvalue A = H(c)/u is genuinely nonlinear. Our next assumption 
is then: 

/" 7^ 0, everywhere 

For Temple systems (see the next section for a definition) it is well known that there is 
no blow up in then we have to investigate if the PSA system is in this class or not: 



8 



this is done in the next section. This will lead us to state this supplementary assumption to 
avoid the Temple class: 



These four hypothesis are not natural, in the sense that they are not generally satisfied, but 
some applications are concerned, as the case of an inert gas associated with ammonia or 
water vapor for instance. 

3 Temple class and PSA system 

It is well known ([5JQ]) that blow up cannot occur for Temple systems with L°° data. The 
aim of this section is thus to investigate the conditions for the PSA system ©-([H]) to be or 
not to be in the Temple class. Of course, in the sequel, we must avoid the Temple class to 
construct our blow up solution. 

About Temple systems, on may consult [2H1 [29j EH [15j [25]. Let us recall the definition 
of the Temple class ([29]). 

Definition 3.1 (Temple class) A system of conservation laws in one space dimension is 
in the Temple class if it is strictly hyperbolic and satisfies the following properties: 

1. there exists a system of coordinates consisting of Riemann invariants, 

2. shock and rarefaction curves coincide. 

For PSA system we have the following criterion: 

Lemma 3.1 PSA system (E|)-(0|) is a Temple system if and only if d x W — for the entropy 
solution of any Riemann problem U6}), (|2ip. 

Proof : as seen above, PSA system is strictly hyperbolic and satisfies the first condition 
of the definition. There is a linearly degenerate eigenvalue (zero) and the rarefaction curves 
are only given by d x W = (see (fH|) or [5]). 

Notice that the Riemann problem for PSA system has a unique piecewise smooth entropy 
solution. If PSA system is in the Temple class, the second condition in the definition ensures 
that d x W = everywhere. The converse is immediate. □ 

Let us highlight a sufficient condition for the PSA system to be in the Temple class: 

Proposition 3.1 If G" = then PSA system is a Temple system. 

Proof : according to (IT51) . if G" = we have d x (uG(c)) = and from Lemma UTTl PSA 
system is a Temple system. □ 

We precise now some simple examples where G" = 0. 
Proposition 3.2 If the two isotherms are linear then G" = 0. 

Proof : we have G = exp(g), G' = g'exp(g), G" = {g" + 5 ' 2 )exp( 5 ). With g' = -ti/H, 
H = 1 + q[ — ch'(c) and h = q\ + qi we get 



PSA system {5p-([6p is not a Temple system 



(HA). 



G" = 



exp(ff) 
H 2 



i-q'i-qZ-qhi+qWO- 
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For two linear isotherms we have q'{ — q 2 = and then G" = 0. □ 
More generally we have: 

Lemma 3.2 Assume (H2), then G = if and only if there exist two real constants a, f3 
such that 

aqi + (a - I) q 2 + c + /? = 0. 

Proof : as seen above, G" — if and only if q" + q 2 + q 2 q[ — q' 2 Qi = 0. Setting 
hi = q[, h 2 = q 2 this condition writes (h\ + h 2 )' — (hi + h 2 )'h 2 ~ h' 2 (hi + h 2 ) or also 

(hi + h 2 )' = (h 1 + h 2 yh 2 -h' 2 (h 1 + h 2 ) 
(hi + h 2 y~ (hi + hrf 

Thus G" = if and only if there exists a constant a such that 

-1 h 2 

= -\- ot, 

hi + h 2 hi + h 2 

i.e. —1 = — q' 2 + oi(q'i + ^2)7 which gives immediately the condition. □ 

In the particular case where a gas is inert, say the first one (qi = 0), we can give a very 
simple criterion: 

Proposition 3.3 Assume qi = 0. Then PSA system {Sp-jOl) is a Temple system if and only 
if the "active" isotherm is linear (q 2 — 0J. 

Proof : if q 2 = and qi = we have G" = thanks to Proposition 13.21 Then the 
system is in the Temple class thanks to Proposition 13. II 

Reciprocally, assume that the system is in the Temple class: following Lemma 13.11 
d x (uG(c)) — holds for the entropy solution of any Riemann problem (|16|) . Let c_, c + , u + 
be the data of the Riemann problem: the unknown is w_ and using the Rankine-Hugoniot 
conditions we have (see [5]): 

u -1Z(c c^u^-i u +7(c-,c+) for a shock 

U - ~ K{C -' C +' U+) ~ { u+ G(c + )/G(c_) for a rarefaction, [Z/) 

where 

[c + qi(c)} - c + [h] a + h- 

7(c-,c+) - 



[c + qi (c)] — c_ [h] a + h 



+ 



If] 

with a = + 1, [ c ] = c + - c-, [/] = /+ - /- = f(c+) - f(c-) (and so on). 
W 

Since the system is in the Temple class, shock and rarefaction curves coincide, thus we 

have: — — — — = — = -7^, then —a [G] = [hG]. We can assume g(0) — 0, then choosing 
a + h + u + G- 

c_ = and setting c = c+, we get [G] = G - 1, [hG] = hG - K with K = h(0) > and 
a = 1 - h(c). Thus, the equality -a [G] = [hG], reads G = 1 + K - h. 

-hi 

Differentiating this relation we get g' G — —h', then G — = H. Differentiating 

9 

H' —h! 
G = H we get g' = — , but in view of (|13l) we have g' = — - thus h' = —H' = h! + ch" and 

finally h" = i.e. q' 2 ' = 0. □ 
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Remark 3.1 For an inert gas (for instance the first one: q\ = 0) and an active gas with 
strictly convex isotherm (q^ > 0), PSA system is not in the Temple class. For example, it is 
the case if the active gas is the ammonia or the water vapor. For other examples, see JH [7]/. 

It is now clear that (H4) can be satisfied in some physically relevant cases. 

4 The blow up solution 

In this section we present and prove our main result: the existence of an entropy blow up 
solution. It is organized as follows: first we state and comment this result, next we give a 
basic step allowing to increase the velocity through two non interacting waves, a shock and 
a rarefaction, then we repeat this process to build a piecewise Lipschitz entropy solution 
with growing velocity. Lastly, a sequence of such solutions leads to a blow up (in a bounded 
space-time domain). 

4.1 Main result 

Our main result needs four assumptions explained in subsection 12.31 

Theorem 4.1 (Blow up for velocity) Assume (HI), (H2), (HS) and (H4). Then for any 
T > 0, Xqo > 0, there exists a set of initial and boundary data and a corresponding weak 
entropy solution on [0, T] x [0,Xoo[ of PSA system such that 

II u IU°°((0,T)x(0,X«,)) = +°°- 

Actually the solution built to prove this theorem has special features. 

• The velocity only blows up at the boundary {t — 0} when x — > X^, with arbitrary 
small data. 

• The system is strictly hyperbolic, but at ( t = 0, x = Xoo) it becomes degenerate 

hyperbolic because A = tends to when u tends to +oo. Then the boundary 

u 

becomes twice characteristic. 

• VX GlO^X^l u,ce L°°(0, T; BV(0, X)) n L°°(0, X; BV(0, T)). The concentration c 
remains bounded while u blows up. 

• Let f2 be a neighborhood of the critical point (t =Q,x= Xoo) such that ft C [0, T] x 
[0,Xoo]. Outside Q, we can prove that (u, c) has a piecewise smooth structure, so 
the blow up occurs only at the boundary. Indeed, there is an accumulation of wave- 
interactions near (t = 0, x = Xoo). 

• To build such solutions, we have necessarily to choose the boundary concentration 
Cq(x) = c(0,x) not in BV^^Xoo). Else there is no blow up, see |U|S]- 

4.2 A mechanism to increase velocity 

In this section, we present the main ingredient for the construction of the blow up solution. 
It consists in the resolution of two consecutive boundary Riemann problems which leads to 
increase the velocity without increasing the concentrations. 
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U C+ 



Figure 4: the two boundary Riemann problems 

For the first boundary Riemann problem the data are (c_, c+, uq) chosen in such a way 
that the solution is a shock wave and we set (see U\ — 7£(c_, c + , u ). For the second 

problem the data are (c+, c_, iti), the solution is necessarily a rarefaction wave and we set 
u 2 = lZ{c + , C_, Ml). 

The sign of G" is crucial but not those of /" and h! in assumptions (H3), (H2). Indeed, 
exchanging the labels of the two gas changes signs of /" and h! . In the sequel, we fix 
these signs: 

/" < 0, ft.' < everywhere 

but it is not restrictive and all variants work. 

With these assumptions we must take c_ < c+ (see Propositions 12.21 and [3]). 
We introduce now the amplification coefficient R defined by 

u 2 = Ruq. 

We are going to show that this coefficient only depends upon (c_, c+). 
Lemma 4.1 The amplification coefficient satisfies R = R(c_,c + ). 
Proof : according to (|2"2")l . for the shock the Rankine-Hugoniot condition gives 

ui = u 7(c_,c+) 

and for the rarefaction, we have d x (uG(c)) = then u 2 G(c+) = u\ GicJ). Finally we have 

U2 = ^ C ~| l{c-,c + )u = R(c_,c+)u - 
G(c+) 

□ 

The following lemma justifies the label "amplification" . 
Lemma 4.2 We assume (HI), then the amplification coefficient satisfies R > 1. 
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Proof : since G" < then d x (uG(c)) > (entropy condition (JTSJ) ) . Through the shock 
wave we have u\ > uq, c_ < c + and uq G{c + ) < U2 G(c + ). Finally, u-i — Huq > uq then 
R>1. □ 

Notice that for the shock curves we have c_ < c+ and u is not monotonic along the 
process because uq < u% < u\. Actually we need a strict amplification i.e. R > 1. We show 
that it is true for almost every choice (c_, c+) with c_ < c+, provided the isotherms are 
analytic as we will see below. We proceed in two steps. 

Lemma 4.3 R(c_, c+) = I if and only if the system is in the Temple class. 

Proof : assume that R(c_, c+) = 1 then U2 = uo- This is true for any (c_, c + ) with c_ < c+ 
and thus the shock curve connecting (c + , uq) to (c_, ui) and the rarefaction curve connecting 
(c_, ui) to (c_|_, U2 = tio) are the same, but traveled in opposite directions. Finally shock and 
rarefaction curves coincide i.e the system is in the Temple class. 

Reciprocally, if the system is a Temple system we have d x (uG(c)) = 0, then 

u a G(c+) = Hi G(c_) = u 2 G{c + ). 

finally uq = Ui and R(c_, c + ) = 1. □ 



Lemma 4.4 We assume (HI), (HS) and (HA). 

If the function R is analytic then R(c_, c+) > 1 almost everywhere in the domain c_ < c+. 

Proof : thanks to (i?l) and lemma [4T^1 we have R > 1. Since R is analytic then either 
R(c_,c + ) = 1 always, either R(c_,c + ) > 1 almost everywhere. With the assumption (HA), 
PSA system is not a Temple system then, thanks to the previous lemma, R > 1 almost 
everywhere. □ 

Notice that the commonly used isotherms are analytic, then R also. With Remark 13.11 
we have some relevant examples where the situation R > 1 occurs. 



4.3 A solution with growing velocity 

Following the mechanism introduced in the preceding subsection we impose constant initial 
data and piecewise constant concentration at the boundary such that Riemann problems at 
the boundary are alternatively solved by a shock or a rarefaction. 

Let N > 1 be a fixed integer, = xq < X\ < ■ ■ ■ < X2N-1 < A = X2N, < c < c < 1 
such that R(c, c) > 1 and u > 0. We solve PSA system ([S)-© with the following data for 
<t <T, < x < X, k = 0,1,--- ,N-1: 

c, 

r -t U(h (23) 

I C if X 2 k < X < X 2 k+1, 
\ C if X 2 k+1 < X < X 2 k+2- 

With these data we prove existence and uniqueness of a piecewise Lipschitz entropy 
solution: 



c(t,0) = 
u(t,0) = 

c(0,x) = 
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Theorem 4.2 We assume (HI), (H2), (H3). With data \23\) . there exists a unique weak 
entropy solution in the class of piecewise Lip schitz functions. Furthermore, this solution has 
exactly N shock curves on [0,T] x [0, A]. 

The uniqueness is an important feature to extend the previous entropy solution until the 
blow up. 

Proof : to construct the solution we use the Front Tracking Algorithm as in [7] (see 
[H] for general references). There exists a domain Z (see Fig. [5]), a neighborhood of the 
boundary {t — 0} where boundary Riemann problems do not interact. In this domain Z we 
have an explicit solution. Let us denote by Uk the value of u(x, + ) when Xk-\ < x < Xk for 
a given k > 0. 

Since /" < and hi < we have N shocks emerging from ((x 2 k,t = 0))k=o and N 
rarefactions from ((x 2 k+i,t = O))^ 1 . Furthermore u 2 k < u 2k +i, u 2 k < u 2 k+ 2 — Hu 2k and 

u 2k = R k u . (24) 

Take 5 > to discretize the rarefactions as in [7] and let shocks, rarefactions and contact 
discontinuities interact to obtain an approximate entropy solution on [0,T] x [0, A]. 

As in a scalar conservation laws with piecewise constant data, no new shock can appear 
but shocks can disappear. It is a consequence of wave- interactions studied in [7]. We recall 
the following results concerning interactions: 

• if a rarefaction interacts with a shock then we have a contact discontinuity and a 
rarefaction or a shock. 

If the shock and the rarefaction have the same strength (i.e. the same absolute variation 
of the concentration through the simple wave) we have only a contact discontinuity; 

• if two shocks interact we obtain a contact discontinuity and a stronger shock; 

• if a shock interacts with a contact discontinuity, we obtain a contact discontinuity and 
a shock with the same strength; 

• if a rarefaction interacts with a contact discontinuity, we obtain a contact discontinuity 
and a rarefaction with the same strength. 

So, we can follow each shock, more precisely: 

• if a shock disappears after an interaction with a stronger rarefaction, by convention 
we follow the characteristic speed associated with the eigenvalue A = H(c)/u, and the 
strength of the shock is 0; 

• if two shocks interact, they become a single shock. But, we consider this shock as two 
shocks to define without ambiguity the parametrization of the N shocks issuing from 
((x 2k ,t = 0))f =1 by N functions sf : [0, X] -> [0,T]. 

The entropy inequality d x u G(c) > yields a lower bound inf u > (see [4]): indeed we have 

uG(c) > u b (t) G(c b (t)), G(c b ) > inf G > 0, < G(c) < supG and u b (t) > a > 0. Then we 

[o.t] [0,1] 

SUPfQ -n H 

get an upper bound for A = H(c)/u: < A < A = — • Thus the speeds of the shock 

mf u 

dsl 

curves — — are bounded by A. 

dx 
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Thus, we have N shock curves uniformly Lipschitz. When 5 goes to 0, Ascoli's Theorem 
gives N Lipschitz curves (sfc)fcLu U P to a subsequence. We denote by U 5 — (c s ,u s ) the 
approximate solution issued from the Front Tracking Algorithm with S > fixed. Extracting 
again a subsequence (6 n ), the sequence U n = U 5 " converges towards an entropy solution 
U = (u, c) as in [U G3 EJ. Let us denote by S the set of shocks emerging from the boundary. 

N 

S = (J T k where T k = {(x,s k {x)), x e [0,X]} 

k=l 

and , S s with a similar definition for the approximate solution II s . Let us show that c has 
no singularity except on S. 

Outside S the concentration c is necessarily continuous: notice that c has BV regularity 
(see [5J [7]), then c is discontinuous if and only if c has shocks. Furthermore (c,u) is an 
entropy solution hence dtc > through the shocks. So it suffices to prove that dtc < 
outside S to conclude. 

For this purpose, take any open set f2 such that S D fi = 0. Using the uniform limit of 
the shock curves, we have no shock for U n in f2 and n sufficiently large: n > n\. So, in Q, 
any jump of c satisfies [c 5n ] < since it is in a rarefaction domain for n > ri\. That is to 
say dtC Sn < in fi for n > n\. Passing to the limit we get dtc < 0. Then, c is continuous in 
f2. This is true for any open set outside S, so c is continuous except on S. 

Since c is continuous outside a finite number of shocks, we can obtain c and u by the 
method of characteristics, as in ;7i . More precisely, it was proved in [7J that any entropy 
solution with a continuous concentration is unique and can be constructed by the method 
of characteristics. Since A is bounded, the characteristics are Lipschitz continuous, so c is 
Lipschitz too in the complementary of S, which prove the piecewise Lipschitz structure of 
the concentration for the entropy solution. Notice that u is less regular than c because there 
are contact discontinuities emerging at each interaction between two waves. 

Uniqueness: the piecewise Lipschitz structure allows us to prove uniqueness. It is a simple 
extension of a theorem proved in [4] for piecewise smooth solutions. For completeness, we 
give the key points of this proof. Under the first shock, namely in the domain Zq (see Fig. 
E]), (u,c) is unique and determined by the method of characteristics. In our case, (it, c) is 
constant in but we do not use this fact to give a general argument and to use it for the 
other domains. 

Over the first shock and under the second one, namely domain Z\, the characteristics are 
well defined if we know u on the upper side of Y\ since d x u G(c) = on Z\. Indeed we have 
the following equations in variables (u, c) for some functions F and S: 

d x W = 0, uiZqUZ-l, 

d t c + W 0(c) d x c = 0, in Z UZ U 

d t x{t,s,x) = W(f>(c), iaZoUZ u 

Wi-Wq = F{cq,c{), onTi, 

— = 5(Wo,c ,ci), on Ti 

with W(t, 0), c(t, 0), c(0, x) are given functions and (j)(c) = l/(H(c)G(c)). 

The trick in 4J is to work in variable t. In the same way, t is an admissible parameter 

for the shock curves T^: — — is well defined since we have BV estimates for lnw, thus A 

dt 
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avoids zero. Moreover W has jumps with respect to t (contact discontinuities for u) and the 

product dtc may not be defined. 

W q>{c) 

The first two equations are equations for the Riemann invariants W and c. The third 
equation is the equation of characteristics associated to the eigenvalue A where x(£, t, x) = x. 
The fourth equation is the Rankine-Hugoniot equation, where Cj and Wi are the traces of 
c in Zi on T\. This relation gives us W in Z\ if we know Ti, cq and c\. Indeed Ti is an 
unknown curve, but the last equation determines it with an ordinary differential equation. 




Figure 5: first interactions and free domain Z 

At t = 0, we know Wo, cq, a so we know Ti and W\ locally near the corner for all t <r, 
and then for all x. We use a local uniqueness theorem as the Cauchy-Lipschitz Theorem. 
Notice that when we know W\ for t < r, the characteristics are also defined for all x thanks 
to d x W — 0. So, we have a uniqueness result in Z U Ti U Z 2 for t < t and we extend the 
uniqueness result until t = T. Next, we use the same technique for the second shock curve 
and so on. 

Now, as in @], U = (c,u) is unique. So, the whole sequence II s tends to U. □ 
4.4 Blow up of the velocity 

Now, we use data (|23l) with N = +oo i.e. (xk)k&i is an increasing sequence such that 
lim Xk = Xoo. Then, the concentration remains bounded at the boundary but the BV 

k— >+oo 

norm of c and the L°° norm of u blow up at the boundary t = 0. Near (t = 0, x = Xoo), there 
is an accumulation of oscillations for the concentration and an accumulation of interactions 
between shocks and rarefactions. 

Let Z be the subset of [0, T] x [0, X^l defined in the previous subsection: it is a neighborhood 
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of the vertical segment {0} x [0, Xoo [. Indeed, the first interaction takes place at (t 1 , x 1 ) : since 
a contact discontinuity propagates vertically, we have first to cut the set {t > t x ,x > x 1 }. 
We do the same work, for all first interactions of the 2N Riemann problems issuing from the 
boundary. 

Then Z = {0 < x < Xoo, < t < z(x)} where the function z : [0,Xoo] -> [0,T] is 
piecewise constant on [0,X] for any X < X^, non increasing, positive before X^, z(0) = T, 
z(X oo ) = 0. 

On Z, we exactly know the solution and 

lim u(t, x) = +oo. 

Proof of Theorem I4.lt let < X < X^. By previous construction, we get a unique 
piecewise smooth entropy solution U on [0,T] X [0, X]. If < X < Y < I M1 we get U Y 
and by uniqueness, U x = U Y on [0, T] x [0,X]. So, as for an ordinary differential equation, 
we can consider the unique maximal solution on [0, T] x [0, X^. Before interaction between 
the solutions of the Riemann problems issuing from the boundary, we know explicitly U in 
Z. From geometric increasing of u(0, x) when x —> X^, we get a blow up for u at t = 0, 

X Xqq . 

Furthermore, the characteristic slope 1/A goes to infinity near {0} x {X^} but outside a 
suitable "triangular neighborhood" of this point we get a determination domain where the 
solution has BV regularity and there is no blow up. This determination domain has the form 

T = {(f, x);Q<t<T, <x < min(X + ^IJ 

A 

for any X in ]0, X tX3 [. □ 



5 Conclusion 

To conclude this paper, let us give some comments about this topic. 

• It is a new physical example of hyperbolic partial differential equations exhibiting a 
blow up. However, this chemical model becomes not physically relevant for too large 
velocity, u » 1. 

• It follows from our blow up example that general L°° theory for PSA system only with 
co,Ch,lnufc e L°° is impossible. But, since our example does not blow up at finite 
distance of the boundary, a natural conjecture is: "if lnu remains bounded at the 
boundary t — 0, there is no blow up". 

• We also conjecture that our solution can be defined for x > X^ as an entropy lo- 
cally piecewise smooth solution outside the boundary, with infinite velocity on {0} x 
[Xoo,+oo[. 

Finally we give some comparisons with other blow up solutions for 2x2 systems. The 
reader can also consult [37] and the memoir |26j . 

• There is a very simple and impressive example in [311 132) , with two Burgers equations 
linearly coupled on two boundaries. This solution blows up in finite time T for all x. 
In contrary, for our example, there is only one boundary and the solution only blows 
up on this boundary. 
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• As in |16] , it is an example exhibiting a blow up of amplitude where at the same time 
the system looses strict hyperbolicity. This is a special feature of 2 x 2 systems, since 
for larger systems, at least 3x3 systems, this is well known that blow up can occur 
without loosing the strict hyperbolicity. 
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